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Abstract 

We present a review of recent advances in the study of many-body effects in mag- 
netic nanoparticles. Considering classical spins on a lattice coupled by the exchange 
interaction in the presence of the bulk and surface anisotropy, we investigate the 
effects of finite size, free boundaries, and surface anisotropy on the average and local 
magnetization for zero and finite temperatures and magnetic fields. Superparam- 
agnetism of magnetic particles necessitates introducing two different, induced and 
intrinsic, magnetizations. We check the validity of the much used relation between 
them within different theoretical models. We show that the competition between 
the exchange and surface anisotropy leads to spin canting dependent on the orien- 
tation of the average magnetization with respect to the crystallographic axes and 
thus to a second-order effective anisotropy of the particle. We have also investigated 
the switching mechanism of the magnetization upon varying the surface anisotropy 
constant. Some cases of more realistic particles are also dealt with. 

Key words: Fine particle systems, surface anisotropy, classical spin models 
PACS: 75.50.Tt - 75.30.Pd - 75.10.Hk 



1 Introduction 



Magnetic nanoparticles, or nanoscale magnetic systems, have generated con- 
tinuous interest since late 1940s as the investigation of their properties turned 
out to be challenging from both scientific and technological point of view. In 
1949, in a pioneering work [1], Neel set the pace towards understanding of 
the magnetic behavior of nanoparticles, leading to an important development 
of fundamental theories of magnetism and modeling of magnetic materials, 
as well as remarkable technological advances, e.g., in the area of information 
storage and data processing, fostering the development of magnetorecording 
media with increasingly higher densities. Nanoparticles, as compared to bulk 
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materials, possess very important novel properties such as enhanced rema- 
nence and giant coercivity, as well as exponentially slow relaxation at low 
temperature due to anisotropy barriers, which ensures great stability of the 
information stored. However, nanoparticles become superparamagnetic [1] at 
finite temperatures for very small sizes, and this is an impairment to the in- 
formation storage. On the other hand, discovery of superparamagnetism that 
results from thermally activated crossing of the anisotropy energy barrier by 
the magnetic moment of the particle opened a rich area to the application of 
nonequilibrium statistical mechanics. There are mainly two types of nanopar- 
ticle samples: i) Assemblies of, e.g., cobalt, nickel or maghemite, nanoparticles 
with volume distribution and randomly oriented axes of magnetocrystalline 
anisotropy; ii) Isolated single particles of cobalt or nickel that can be probed 
by the micro-SQUID technique [2]. While most of the experiments have been 
done on assemblies, isolated particles are more important both as units of 
information storage and as a physical system. 

In the investigation of the static properties of magnetic nanoparticles a great 
deal of work up to date has been based on the Monte Carlo (MC) technique. In 
addition to numerous simulations of the Ising model, this technique has been 
used with the more adequate classical Heisenberg model to simulate idealized 
isotropic models with simple cubic (sc) lattice and spherical shape in Ref. [3]. 
Magnetic nanoparticles with realistic lattice structure were recently simulated 
in Ref. [4] taking into account the surface anisotropy (SA) and DDI. 

On the other hand, in most of theoretical approaches to the dynamics of a 
small magnetic particle the latter is considered as a single magnetic moment. 
This is the one-spin approximation that is only valid for particles that are 
not too large and thus are single-domain, and not too small to be free from 
surface effects. Letting apart spin tunneling (see, e.g., Ref. [5]) that becomes 
important for extremely small sizes such as those of molecular magnets, the 
magnetic moment can overcome the anisotropy-energy barrier and thus reverse 
its direction, at least in two ways 1 : Either under applied magnetic field 
that suppresses the barrier, or via thermal fluctuations. The former, at zero 
temperature, is well described for particles with the uniaxial anisotropy by 
the Stoner-Wohlfarth model [7]. Thermally activated crossing of the energy 
barrier is described by the Neel-Brown model [8] and its extensions (see, e.g., 
Refs. [9]). At elevated temperatures, rotation of the magnetization in materials 
with strong anisotropy is always accompanied by changing its magnitude. This 
results in a shrinking of the Stoner-Wohlfarth astroid as described by the 
modified Landau theory [10], and (qualitatively) confirmed by experiments 
[2]- 



1 It has quite recently been shown [6], experimentally and theoretically, that effi- 
cient magnetization switching can be triggered by transverse field pulses of a dura- 
tion that is half the precession period. 
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Both Stoner-Wohlfarth and Neel-Brown models have been confirmed by ex- 
periments on individual cobalt particles [2]. However for magnetic particles 
with strong surface anisotropy magnetization switching occurs as the result 
of successive switching of individual (or clusters of) spins inside the particle 
[11]. Such deviations from the one-spin approximation have been observed in 
metallic particles [14], [16], and ferrite particles [12], [13]. Deviations from 
the one-spin approximation and temperature effects lead to the absence of 
magnetization saturation at high fields [14], [15], [16], shifted hysteresis loops 
after cooling in field, and field dependence of the magnetization at very low 
temperatures. The latter effect has been clearly identified in dilute assemblies 
of maghemite particles [17] of 4 nm in diameter. In addition, aging effects 
have been observed in single particles of cobalt and have been attributed to 
the oxidation of the sample surface into antiferromagnetic CoO (see [2] and 
references therein.). It was argued that the magnetization reversal of a ferro- 
magnetic particle with antiferromagnetic shell is governed by two mechanisms 
that are supposed to result from the spin frustration at the core-shell interface 
of the particle. Some of the above-mentioned novel features are most likely due 
to magnetic disorder at the surface which induces a canting of spins inside the 
particle, or in other words, an inhomogeneous magnetic state. This effect was 
first observed with the help of Mossbauer spectroscopy by Coey [18] and later 
by Morrish and Haneda (see [19] for a review), and later by Prene et al. [20] 
(see also the recent article [21]). To sum up, the picture of a single-domain 
magnetic particle with all spins pointing into the same direction is no longer 
valid when one considers the effect of misaligned spins on the surface, which 
makes up to 50% of the total volume in a particle of 4 nm in diameter. 

One of our goals is thus is to understand the effect of surfaces on the thermo- 
dynamics and magnetization profiles in small systems, and subsequently upon 
their dynamics. This requires a microscopic approach to account for the local 
environment inside the particle, microscopic interactions such as spin-spin ex- 
change, DDI, and the magneto-crystalline bulk and surface anisotropy. As this 
task is difficult owing to the large number of degrees of freedom involved, one 
has to gain a sufficient understanding of static properties before proceeding 
to the dynamics. There are three main effects that distinguish magnetic par- 
ticles from bulk magnets and that were investigated in a series of our recent 
publications: 

(1) Finite-size effect in isotropic magnetic particles with idealized periodic 
boundary conditions. The spin-wave spectrum of such particles is discrete 
and there is the mode with k = that corresponds to the global rotation 
of the particle. As a result, the standard spin-wave theory fails and one 
has to distinguish between induced and intrinsic magnetizations. 

(2) Boundary effect, i.e., pure effect of the free boundary conditions at the 
surface in the absence of the surface anisotropy. This effect leads to the 
decrease of the particle's intrinsic magnetization and it makes the latter 
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inhomogeneous at T 7^ 0. 
(3) Effect of surface anisotropy that changes the ground state of the particle 
and makes the intrinsic magnetization inhomogeneous even at T = 0. 

This paper is organized as follows: In Sec. 2 we describe the Hamiltonian and 
introduce the basic notions of the induced and the intrinsic magnetizations. 
In Sec. 3, we first study inhomogeneities in small magnetic particles of a box 
shape induced by pure boundary effects in the absence of surface anisotropy, 
i.e., the effect of free boundary conditions (fbc), and compare their influence 
with that of the finite-size effects in an idealized model with periodic boundary 
conditions (pbc). We consider spins as D-component classical vectors. At first 
we present analytical and numerical results in the whole range of tempera- 
tures in the limit D — > 00 where the problem simplifies. Then for the classical 
Heisenberg model, D — 3, at low temperatures we formulate the modified 
spin-wave theory accounting for the global rotation of the particle's magneti- 
zation. The results are compared with those of our MC simulations. In Sec. 4 
we consider round-shaped systems and include surface anisotropy. In the case 
of the SA much weaker than the exchange interaction we study the problem 
perturbatively in small deviations from the perfectly ordered, collinear state. 
Then we investigate the hysteretic properties and the behavior of the mag- 
netization as a function of temperature and applied field by MC simulations. 
The last section summarizes the results and points out open problems. 



2 Basic relations 

2.1 The Hamiltonian 

Within the classical approximation it is convenient to represent the atomic 
spin as the three-component spin vector Sj of unit length on the lattice site 
i. We will consider the Hamiltonian that in general includes the exchange 
interaction, magneto-crystalline anisotropic energy, Zeeman energy, and the 
energy of dipolar interactions (DDI) 



where no = g^sS and S is the value of the atomic spin. For materials with 
uniaxial anisotropy 7i an in Eq. (1) reads 




ij 



/UqH • Sj + Han + HbBI, 



(1) 



ftin ni) = -£^(s;.e0 2 , 



(2) 
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with easy axis e« and constant Ki > 0. This anisotropy model can be used to 
describe the surface effect if one attributes the same easy axis and the same 
anisotropy constant K c for all core spins and different easy axes and anisotropy 
constants to surface spins. Within the simplest transverse surface anisotropy 
(TSA) model all surface spins have the same anisotropy constant K S) whereas 
their easy axes are perpendicular to the surface, see. e.g., Refs. [22,23,24] . 
More realistic is Neel's surface anisotropy (NSA) model [25] , 

z« 

^an ^ = — 5Z( S * ' e *i) ' e *i = r *j/ r U> r *i = r i ~ r j> (3) 

i j=l 

where Zi is the coordination number of site i that for the surface atoms is 
smaller than the bulk value z and is the unit vector connecting the site % 
to its nearest neighbors j. One can check that for the simple cubic (sc) lattice 
the contributions from the bulk spins in (3) ~ sf = 1 are irrelevant constants. 

For materials with magneto-crystalline cubic anisotropy 7i an reads 

w£ = -k<* e (4 + 4 + 4) • (4) 



For > 0, the energy TiS^ has minima for six orientations of the type [100] 
and maxima for eight orientations of type [111]. 

If one discards TiS^} and 7Yddi! the Hamiltonian (1) can be generalized for 
.D-component spin vectors. This is useful as in the limit D — > oo the problem 
simplifies while retaining important physics, see below. 



2.2 Magnetization of finite systems 

Magnetic particles of finite size do not show magnetic ordering at nonzero 
temperatures at H = as the global magnetization of the particle can assume 
all possible directions (superparamagnetism). It is thus convenient to define 
two magnetizations, m and M, the first being the magnetization induced by 
the magnetic field and the second being a measure of the short-range order in 
the particle. Omitting the factor /xq that can be restored later, we first define 
the magnetization of a microscopic spin configuration 
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where M is the number of magnetic atoms in the system. The thermodynamic 
average of M yields what we call the induced magnetization 



m=(M> = ^£(s,>. (6) 



The intrinsic magnetization is related to the spin correlation function: 
M = yj (M 2 ) 





(?) 



If the temperature is low and there is no surface anisotropy, all spins in the 
particle are bound together by the exchange interaction and M behaves as a 
rigid "giant spin", |M| = M=l. If a magnetic field H is applied, M exhibits an 
average in the direction of H, which leads to a nonzero value of the induced 
magnetization m. For isotropic /^-component vector models (including the 
Ising model, D — 1) [26] the latter is given by the well known formula 

m = MB D (Mx), x = AfH/T, (8) 



where Bjj(x) is the Langevin function [B 3 (x) = cothrr — 1/x for the isotropic 
Heisenberg model and Bi(x) = tanhx for the Ising model]. An important 
question is whether Eq. (8) remains valid at elevated temperatures where 
M = M(T,H). We have shown that the superparamagnetic relation, Eq. (8), 
becomes exact for D — > oo but otherwise it contradicts the exact relation 

. ^ 9 9 dm (D - l)m 

M =m+ — h — . (9) 

dx x 



One also can introduce the local intrinsic magnetization Mj according to 




1^M, = M. (10) 



This quantity is smaller near the boundaries of the particle than in the core 
for the model with fbc because of boundary effects at T > 0. 
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3 Nanoparticle as a multi-spin system: finite-size vs boundary ef- 
fects 

Finite-size magnetic systems with free boundary conditions (fbc) present a 
spatially inhomogeneous many-body problem. In this section we shall only 
deal with finite-size versus boundary effects, leaving the more profound effects 
of the surface anisotropy for the next section. One of the interesting problems 
here is the interplay between boundary effects due to fbc and the "pure" finite- 
size effects. In systems of hypercubic shape, the latter can be singled out by 
using artificial periodic boundary conditions (pbc). The standard mean-field 
approximation (MFA) and spin-wave theory (SWT) are inappropriate for finite 
systems because of the Goldstone mode associated with the global rotation 
of the magnetic moment in zero field. Appropriate improvements include the 
so-called D — > oo model (see Sec. 3.1) operating at all temperatures and 
the modified spin-wave theory for finite magnets at low temperatures (see 
Sec. 3.2). Also the MC routine should incorporate global rotations of spins, in 
addition to the Metropolis algorithm of individual spin rotations. 



3.1 D — > oo model 
3.1.1 The model 

One can improve upon the MFA by taking into account correlations in a wide 
temperature range for bulk and finite magnets by replacing 3-component spin 
vectors in (1) by /J-component ones and taking the limit D — > oo. This model 
was introduced by Stanley [26] who showed that in the bulk its partition 
function coincides with that of the exactly solvable spherical model (SM) [27]. 
On the other hand, for spatially inhomogeneous and anisotropic systems the 
D — > oo model is the only physically acceptable model of both (see, e.g., Refs. 
[28]). So far, the D — > oo model was only applied to spatially inhomogeneous 
systems in the plane geometry [28]. In Ref. [10] we extended it to finite box- 
shaped magnetic systems with free and periodic boundary conditions. In the 
MFA the Curie temperature of the D-component model is T C MFA = Jq/D, 
where Jo is the zero Fourier component of J^. It is convenient to use T C MFA as 
the energy scale and introduce the dimensionless variables 

9 = T/T^ FA , h = H/J , \ l3 = J l3 /J Q . (11) 

For the nearest-neighbor (nn) interaction Jij with z neighbors, \j is equal 
to 1/z if sites i and j are nearest neighbors and zero otherwise. In the bulk 
the D — > oo model is described by two coupled nonlinear equations for the 
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magnetization m and the so called gap parameter G: 

m = _^_, m > + eGP(G) = l, P( G)= /^ 1 (12) 



where P{G) is the lattice Green function and Ak = (cos k x + cos k y + cos k z ) /3 
is the Fourier transform of The Curie temperature is defined by G = 1 
and is T c = T™ FA /W, i.e., # c = 1/W, where W = P(0) is the Watson integral 
(W = 1.51639 for the sc lattice). The system of equations describing the 
inhomogeneous D — > oo model can be obtained using the diagram technique 
for classical spin systems [29,30] in the limit D — > oo and generalising the 
results of Ref. [28] for spatially inhomogeneous systems to include the magnetic 
field h = he z . This is a system of equations for the average magnetization 
rrii = (s Z i) directed along the field, gap parameter Gj, and correlation functions 
for the remaining "transverse" spin components labeled by a > 2, i.e., Sij = 
D (s ai s a j) (all transverse correlation functions are the same). This system of 
equations has the form 

^VijUlj = h, J2 V v s ji = s a + m i = !> ( 13 ) 



where X>y = G^ l 5ij — A^- is the Dyson matrix and 5^ is the Kronecker symbol. 
Solving this system of equations consists in determining rrii and s^- as functions 
of Gi from the first two linear equations and inserting the solutions into the 
third nonlinear equation (the constraint equation) that leads to a system of 
nonlinear equations for all Gj that is in general subject to numerical solution. 
Combining Eqs. (13) results in m 2 + 6m/ (Afh) — M 2 = that yields 

m = M 2 Mh/6 = = mb jk[ MH i t \ ( 14 ) 

1 + ^/1+ (2MMh/6) 2 



where 



M 



m2 + TpE%' B °°(0 

ij 



l + y/l + {2i/Df 



(15) 



and is -Boo(0 the Langevin function for D 3> 1. Alternatively Eq. (14) can be 
derived from Eq. (9) replacing D — 1 =^ D and neglecting dm/dx in the limit 
.D — > oo and then solving the resulting algebraic equation for m. One can 
find in the literature formulae of the type m = M s B(AfM s H/T), where the 
saturation magnetization M s is usually associated with the bulk magnetization 
at a given temperature (see, e.g., Refs. [31]). In our case, Eq. (14) is exact and 
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Fig. 1. Left: Temperature dependence of the intrinsic magnetization M, Eq. (15), 
and local magnetizations Mj, Eq. (10 ) in zero field. Right: Long-range magnetiza- 
tion profile in the direction from the center of the cube to the center of a face at 
temperature 9 = T/T C MFA = 0.189714. 

M = M(T,H) is defined by Eq. (15). For large sizes A/", Eq. (14) describes 
two distinct field ranges separated at H ~ H v where 



v 



TD 
ATM' 



(16) 



In the range H < Hy the total magnetic moment of the system is disoriented 
by thermal fluctuations, m < M. In the range H > H v , the total magnetic 
moment is oriented by the field, m approaches M, and both further increase 
with the field towards saturation (m = M = 1) due to the suppression of spin 
waves in the system. This scenario is inherent to all 0(D) models [31]. 

Having established the superparamagnetic relation, Eq. (14) we are left with 
the problem of calculating M(T,H). For the pbc the solution becomes ho- 
mogeneous and one obtains (12) where in P(G) the integral is replaced by a 
sum over discrete wave vectors [10] , whereas M = ^Jm 2 + 6G/ [Af(l — G)]. 
For the model with fbc analytical solution is only possible at low and high 
temperatures. At 9 1 small deviations from the collinear state with M = 1 
can be described by the modified SWT for arbitrary D, see Sec. 3.2. 



3.1.2 Numerical results and discussion 

The method for solving the D — > oo model consists in obtaining the correla- 
tion functions s^- and magnetization m, from the first two linear equations in 
Eq. (13), substituting them into the third equation of Eq. (13), and solving the 
resulting system of nonlinear equations for the gap parameter numerically. 
Fig. 1 (left) shows the temperature dependence of the intrinsic magnetization 
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Fig. 2. Field dependence of the intrinsic magnetization M and induced magnetiza- 
tion 77i for hypercubic lattices with fbc in three and two dimensions. Dashed line 
is a plot of Eq. (14) in which M(H, T) is replaced by its zero-field value. Bulk 
magnetization in two and three dimensions is shown by solid lines. 

M, Eq. (15), and local magnetizations Mi, Eq. (10), of the 14 3 cubic system 
with free and periodic boundary conditions in zero field. For periodic boundary 
conditions, M exceeds the bulk magnetization at all temperatures. In particu- 
lar, at low temperatures this is in accord with the positive sign of the finite-size 
correction to the magnetization, Eqs. (24) and (26). Local magnetizations at 
the center of the faces and edges and those at the corners decrease with tem- 
perature much faster than the magnetization at the center. One can see that 
below the bulk critical temperature M is smaller than the bulk magnetization. 
This means that the boundary effects suppressing M are stronger than the 
finite-size effects that increase M. This is also seen from the low-temperature 
expression of M given in Eq. (24), see also Fig. 3. Fig. 1 (right) shows the 
magnetization profile in the direction from the center of the cube to the center 
of a face. It is seen that perturbations due to the free boundaries extend deep 
into the particle, whereas the MFA predicts, on the contrary, a fast approach 
to a constant magnetization when moving away from the boundary [3]. This 
is a consequence of the Goldstone mode which renders the correlation length 
of an isotropic bulk magnet infinite below T c . MC simulations of the classical 
Heisenberg model [3], [4] yield a similar result (see Fig. 14 right). 

We have shown that the critical indices for the magnetization at the faces, 
edges, and corners are higher than the bulk critical index (3 = 1/2 for the 
present D = oo model. The critical index at the face (3x is the most studied 
surface critical index (see, for a review, Refs. [32]). The exact solution of Bray 
and Moore [33] for the correlation functions at criticality in the D = oo model 
and application of the scaling arguments yield the value (3\ = 1 (see Table II 
in [32]). Exact values of the edge and corner magnetization indices, (3 2 and j3 3 , 
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seem to be unknown for D = oo. Cardy [34] used the first-order e-expansion 
to obtain /^(a) for the edge with an arbitrary angle a. For a = 7r/2 and 
D = oo in three dimensions the result for the edge critical exponent reads = 
13/8 + 0(e 2 ) = 1.625 + 0(e 2 ). To estimate the magnetization critical indices 
in our model we have performed a finite-size-scaling analysis (see, for a review, 
Ref. [35]) assuming the scaling form M = N~PI U Fm{tN 1 I u ) and plotting the 
magnetization times N 13 ^ vs tN 1 / u . Here v — 1 is the critical index for the 
correlation length in the bulk and r = T/T c — 1, where T c = T^ FA /W is 
the bulk Curie temperature. Our results for the system with N = 10 and 
N = 14 merge into single "master curves" for p 1 = 0.86, (3 2 = 1.33, and 
= 1.79, which have been obtained by fitting M oc N" 13 ^ at T = T c , i.e., 
= 6 C = 1/W. Note that our value 0.86 for the surface magnetization critical 
index fli is substantially lower than the value f3\ = 1 following from scaling 
arguments. This disagreement is probably due to corrections to scaling which 
could be pronounced for our small linear sizes N — 10 and 14. A more efficient 
way for obtaining an accurate value of f3\ is to perform a similar analysis for 
the semi-infinite model. The latter was considered analytically and numerically 
for T > T c and H = in [36]. We also mention the Monte Carlo simulations 
of the Ising model [37] which yield fa = 0.80, (3 2 = 1.28, and (3 3 = 1.77. 

The field dependence of M and m at fixed temperature, as obtained from 
the numerical solution of Eqs. (13) for cubic and square systems is shown in 
Fig. 2. Naturally the numerical results for m confirm Eq. (14) which describes 
both the effect of orientation of the system's magnetization by the field and 
the increase of M in field. Using the zero-field value of M in Eq. (14) leads 
to a poor result for m as shown by the dashed curve for the 10 3 system in 
Fig. 2. The field dependence of the particle's magnetization similar to that 
shown in Fig. 2 for the cubic system was experimentally obtained for ultrafine 
cobalt particles in [16], as well as in a number of previous experiments. The 
curves for the square system in Fig. 2 illustrate the fact that in two dimensions 
thermal fluctuations are much stronger than in 3d, which leads to lower values 
of both M and m at a given temperature. The bulk magnetization in 
two dimensions vanishes at zero field and it thus goes below the intrinsic 
magnetization M in the low-field region. 



3.2 Modified spin-wave theory: Low-temperature properties and the super- 
paramagnetic relation 

As briefly discussed in Sec. 2, and in more details in the previous section, it 
is important to investigate the precise relation between the induced magne- 
tization m of Eq. (6) and intrinsic magnetization M of Eq. (7) for the more 
realistic Heisenberg model, D = 3. To do so, we have developed a finite-size 
spin-wave theory that yields analytical results at low temperatures. We also 
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performed simulations with the improved Monte Carlo technique [38]. 

3.2.1 Spin-wave theory for finite- size magnetic particles 

In the absence of SA, at low temperatures all spins in the particle are strongly 
correlated and form a "giant spin" M defined in Eq. (5) which behaves su- 
perparamagnetically. In addition, there are internal spin-wave excitations in 
the particle that are responsible for M(H, T) < 1 at T > 0. These excitations 
can be described perturbatively in small deviations of individual spins Sj from 
the global direction of M. Thus we use M =Ain with |n| = 1 and insert an 
additional integration over d~M. = Ai D ~ 1 dAidn in the partition function, 

Z = J M^dMdn n dSiS (m-^ £ e~ n/T , (17) 

and first integrate over the magnetization magnitude Ai. Thus we reexpress 
the vector argument of the 5-function in the coordinate system specified by 
the direction of the central spin n : 

6 i M ~Jf^ s ') =6 { M ~Jf ^ n ' s * } ) 6 (a7 ? h ~ n(n ' Sl)] ) - (18) 

Then after integration over M. one obtains 

z= jdnz n , z n = yn^^E^- n ( n -^)]) e_Weff/T ' ( i9 ) 

where Z n is the partition function for the fixed direction n and 
H eS = -(n • H) ^(n • Si) - \ E JijSi • Sj - (D - l)Tln 

i ij 

In Eq. (19), the 5-function says that the sum of all spins does not have a 
component perpendicular to M. This will lead to the absence of the k = 
component of the transverse spin fluctuations. That is, the global-rotation 
Goldstone mode that is troublesome in the standard spin- wave theory for finite 
systems, has been transformed into the integration over the global direction 
n in Eq. (19). Z n is computed at low temperature by expanding 7i ca up to 
bilinear terms in the transverse spin components Ilj = Sj — n(n ■ Sj): 

H eS = E - Afn ■ H + \ E A iJ n i ■ n i> 



.(20) 
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Aij = 



&ij Jij i 



(21) 



where E = — (1/2) J2ij Jij is the zero-field ground-state energy. Next, upon 
computing the resulting Gaussian integrals over II", one obtains 



exp 



-£o+.Afn.H\^_ 1 



T 



(2ttT) 



JV-l 



n 



, 'Ah 



(D-l)/2 



(22) 



where for particles of cubic shape 

A k = A + J k -J , A = J2Aj = (D-l)T/tf + n-H. (23) 



The prime on the product in (22) means omitting the k = mode. Results for 
particles of arbitrary shape can be found in Ref. [38]. In Eq. (22) the expo- 
nential factor corresponds to rigid spins whereas the second factor describes 
spin-wave corrections. The latter makes the angular dependence of Z n more 
complicated. Differentiating Z with respect to H yields the induced magneti- 
zation m, then the intrinsic magnetization M can be obtained from the exact 
relation (9), and the validity of the superparamagnetic relation (8) can be 
checked. 



3.2.2 Induced and intrinsic magnetizations, and superparamagnetic relation 



Now we consider particles of cubic shape (N = N 3 ) at low temperatures. For 
both pbc and fbc, we have obtained the following correction to the magneti- 
zation in zero field [10,38] 



D - 1 W N T D - 1 
M = 1 - t, t= — = -^r-W N 6, 



Jo 



2D 



W, 



N 



where W N is the sum without the k = term. The results for pbc and fbc 
differ only by the values of the discrete wave vectors in Eq. (24) [10]: 



k a — < 



2rm a /N, pbc 
irn a /N, fbc 



0,l,...,iV-l 



(25) 



where a — x,y, z. This subtle difference is responsible for much stronger ther- 
mal fluctuations in the fbc model due to boundary effects. The limit iV — > oo 
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Fig. 3. Lattice sums Wn for cubic systems with free and periodic boundary condi- 
tions. W = 1.51639 is the bulk value for the sc lattice. 

of Wn is the so-called Watson's integral W = -P(O) of Eq. (12). The difference 
between W N and W has different signs for pbc and fbc models [10] 



W N - W 
A N = — ^— = < 



W 



0.90 

9ffi(1.17iV) 

2-kWN 



pbc 
, fbc. 



(26) 



Therefore, Fig. 3 shows that the coefficient in the linear-^ term in Eq. (24) is 
smaller than in the bulk for the pbc system and greater for the fbc system. 
That is, boundary effects suppress the intrinsic magnetization at low temper- 
atures while finite-size effects lead to its increase. Fig. 3 shows that boundary 
effects render a larger contribution than the finite-size effects, making the net 
magnetization well below that of the bulk. 

The field dependence of m and M is defined by the expansion of the lattice 
Green function for small gaps, 1 — G <C 1 



Pn{G) 



1 

Jr 



v i 



i 



G\, 



W, 



N 



c N N(l-G), N 2 (l-G) < 1 
c 0V / r^G ? , iV 2 (l-G)>l, 



(27) 



where for the sc lattice cq = (2/7r)(3/2) 3//2 and the numerical results for cjy 
can, for N >> 1, be fitted as 



c N = < 



0.384 - 1.05/iV, pbc 
1.90 - 1.17/JV, fbc. 



(28) 
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The spin-wave gap 1 — G depends on the temperature and field; For H <C Hy 
[see Eq. (16)] the gap approaches its zero-field value while for H ^> Hy one 
has 1 — G = h = Hj J . Thus Eq. (27) defines one more crossover in field, the 
crossover between its first and second lines at 

J J TP 9J 

Hs ~ jpm = m » Hv = mm = mm- (29) 



In the field range H ^> Hs the discreteness of the lattice can be neglected 
and the bulk result AM ~ \fH is reproduced. In the most interesting region 
H <C Hs one obtains 

M=l-t + 2axB(x), a = {D ~^ CN (£)' . (30) 



On the other hand, Eq. (30) describes a crossover from the quadratic field 
dependence of M at low field, x <C 1, to the linear dependence at x >> 1. Note 
that for x 1, where m = M and a rigid magnetic moment would saturate, 
m continues to increase linearly as m = 1 — t + lax. This is due to the field 
dependence of the intrinsic magnetization M. At higher fields there is another 
crossover to the standard spin-wave theory expression for M. Approximate 
expressions for M in the different field ranges are 



f D- 1 (HN T 



M = 1 — t + < 



-c N 



, H<^H V 



2D 

»~\, XHT 
D-l T /H\ 1/2 

— Co T 



(31) 



{-) ,H S ^H « J . 



A simple analysis shows [38] that the superparamagnetic relation (8) is a very 
good approximation for not too small systems, Af ^> 1 in the whole range 
below T c . The deviation from Eq. ( 8) is controlled by the small parameter a 
of Eq. ( 30). Above T c , however, deviations from Eq. (8) are large, except for 
the model with D — > oo. In the close vicinity of T c , there is a crossover to the 
high-temperature form of Eq. (8) given by the function B^x) of Eq. ( 15). 

The modified SWT developed above can be applied to study inhomogeneities 
in the fbc model. The local intrinsic magnetization defined by Eq. (10) shows 
stronger temperature dependence near the boundary than the averaged M = 
1 — t of Eq. (24). The biggest effect of the surface is naturally attained at the 
corners of the cube where Mj = 1 — 8t, at H = [38]. 
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Fig. 4. Comparison of the theoretical and MC results for the field dependences of 
the magnetizations M and m for the Heisenberg model at T = T c /4. 

3.2.3 MC simulations 

Here we apply our Monte Carlo technique that accounts for the global-rotation 
Goldstone mode to compute the induced and intrinsic magnetizations, and to 
investigate the superparamagnetic relation between them, for the Heisenberg 
model with fbc. Our results for Ising model can be found in [38]. 



First of all, in Fig. 4 we compare theoretical predictions of our analytical 
calculations within spin-wave theory for the Heisenberg model with our MC 
results at T = T c /4, where T c = 0.722T C MFA is the actual bulk Curie tem- 
perature. For the small size Af = 5 3 the square-root field dependence of the 
magnetization (third line of Eq. (31)) does not arise and finite-size corrections 
are very important. For M one should use Eq. (30) with t given by Eq. (24) 
with numerically computed Wn = 1-99 and cn = 1-66 for the fbc model [cf. 
Eqs. (26) and (28)]. This yields * ~ 0.119 and a ~ 1.20 x 10" 4 . The theoreti- 
cal dependence M(H) is practically a straight line which goes slightly above 
the MC points. This small discrepancy can be explained by the fact that the 
applicability criterion for our analytical method, t <C 1, is not fully satisfied 
at T = T c /4, and a better agreement is achieved at lower temperatures. For 
comparison we also plot the theoretical M(H) for the model with periodic 
boundary conditions. Here one has = 1.25 and = 0.20, thus t ~ 0.075 
and a ~ 1.45 x 10~ 5 , so M(H) goes noticeably higher and with a much smaller 
slope. The quadratic field dependence of M in the region x < 1 is not seen at 
this low temperature since the value of a is very small and thus much more 
accurate MC simulations would be needed. We also plot in Fig. 4 the analyt- 
ical result for the field dependence of m [38] which favorably compares with 
our MC data. Fig. 5 (left) shows the intrinsic magnetization M and induced 
magnetization m versus the scaled field x = AfH/T for different temperatures. 



16 




Fig. 5. Left: Field dependence of the intrinsic magnetization M and the induced 
magnetization m of the Heisenberg model on the sc lattice with fbc for different tem- 
peratures. Right: Scaled graph for m/M. Theoretical curves B^(x) = cothx — 1/x 
for T <^T C and B OQ (x) for T » T c are shown by solid lines. 

We see that the particle's magnetic moment is aligned and thus m ~ M for 
x > 1, if T <C T c . At T >> T c the field aligns individual spins, and this re- 
quires H >T, i.e., x > A/". The quadratic dependence of M(H) at small fields 
manifests itself strongly at elevated temperatures. 

The results of Fig. 5 (right) show that the superparamagnetic relation of 
Eq. (8) with M = M(T, H) is a very good approximation everywhere below 
T c , for the Heisenberg model. This also holds for the Ising model [38]. On the 
other hand, above T c Eq. (8) with the function B^x) of Eq. (15) is obeyed. 
The difference between these limiting expressions decreases with increasing 
number D of spin components and disappears in the spherical limit (D — > oo). 



4 A nanoparticle as a multi-spin system: Effect of surface anisotropy 

Surface anisotropy causes large deviations from the bulk behavior that are 
much stronger than just the effect of free boundaries. Also SA exerts in- 
fluence upon the coercive field. Here we first consider magnetic structures 
and hysteresis loops at zero temperature induced by a strong transverse sur- 
face anisotropy (TSA). Then for the more physically plausible Neel's surface 
anisotropy (NSA) we investigate analytically and numerically its contribution 
to the effective anisotropic energy of the particle in the case when the NSA 
is weak in comparison to the exchange. After that we study the effect of the 
NSA at finite temperatures using the MC technique. In the end of this section 
we investigate the thermal and spatial dependence of the magnetization of a 
maghemite particle. 
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Fig. 6. Hysteresis for ip = 7r/4, j = 10 2 , N = 10 (M = 360) and different k s . 
4-1 Magnetic structure and hysteresis atT = 0: Transverse surface anisotropy 

Here we study the effect of TSA on the hysteresis loop and the angular de- 
pendence of the switching field. We construct an effective Stoner-Wohlfarth 
(SW) astroid for a single-domain spherical particle with free surfaces, a simple 
cubic (sc) crystal structure, ferromagnetic exchange J, uniaxial anisotropy in 
the core K c , and the TSA of strength K s . Using the Hamiltonian (1) without 
the DDI at T = 0, we solve the coupled Landau-Lifshitz equations (LLE) 
for each spin in the particle [11] until a stationary state is attained. In [39] 
the same method was used for studying hysteresis loops in nanoparticles with 
a random bulk or surface-only anisotropy. Here we address the question of 
whether one can still use the simple SW model for a nanoparticle endowed 
with strong surface effects. We show that it is so as long as K s < J. Otherwise 
switching of the particle's magnetization occurs via the reversal of clusters of 
spins, invalidating the simple SW model. 

We consider K s and exchange coupling on the surface as free parameters since 
there are so far no definite experimental estimations of them, whereas the core 
parameters are taken as for the bulk system. In the sequel, we use the reduced 
parameters, j = J/K c , k s = K s /K c . 

Fig. 6 shows that when k s becomes comparable with j, the hysteresis loop 
exhibits multiple jumps, which can be attributed to the switching of different 
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spin clusters containing surface spins whose easy axes make the same angle 
with the field direction. The hysteresis loop is characterized by two field val- 
ues: One that marks the limit of metastability, called the critical field or the 
saturation field, and the other that marks the magnetization switching, and 
is called the switching field or the coercive field. This progressive switching of 
spins is illustrated in Fig. 7 for simplicity for the non-interacting case. 

h = -1.47 h = 0.47 h = 0.52 
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Fig. 7. Magnetic structures in the middle plane of the particle with N = 360 in the 
TSA model for j = 0, k s = 1, ip = 7r/4 and different fields. 



For small values of k s /j ~ 0.01 the TSA model renders hysteresis loops and 
limit-of-metastability curves that scale with the SW results for all values of the 
angle tp between the core easy axis and the applied field, the scaling constant 
being M c /M < 1, see in Fig. 8. 

For larger values of k s /j, but k s /j < 0.2, we still have the same kind of scaling 
but the corresponding constant now depends on the angle ip between the core 
easy axis and field direction. This is reflected by a deformation of the limit- 
of-metastability curve (see Fig. 9a). More precisely, the latter is depressed in 
the core easy direction and enhanced in the perpendicular direction. However, 
there is still only one jump in the hysteresis loop implying that the magnetiza- 
tion reversal can be considered as uniform. For k s /j > 1, there appear multiple 
steps in the hysteresis loop associated with the switching of spin clusters. It 
makes the hysteresis loop both qualitatively and quantitatively different from 
those of the SW model, as the magnetization reversal can no longer be con- 
sidered as uniform. In addition, in the present case, there are two more new 
features: the values of the switching field are much higher than in SW model, 
and more importantly, its behavior as a function of the particle's size is oppo- 
site to that of the previous cases (compare Figs. 8a and 9b). More precisely, in 
this case one finds that this field increases when the particle's size is lowered. 
This is in agreement with the experimental observations in nanoparticles (see, 
e.g., [14] for cobalt particles). The whole situation is summarized in Fig. 10 
where we plot the critical field h c as a function of k s = k s /j for different values 
of the surface-to- core ratio of the exchange coupling. For large values of k s , 
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Fig. 8. a) (in diamonds) Switching field for k s = 1, j = 10 2 versus the particle's 
diameter N. (in circles) SW switching field multiplied by N c /J\f. b) Astroid for 
k s = 1, j = 10 2 for different values of the surface-to- volume ratio N st = N s /N. 




Fig. 9. a) Astroid for j = 10 2 , N = 360 and different values of surface anisotropy 
constant k s . The full dark line is the SW astroid scaled with N c /M, but the dotted 
line is only a guide for the eye. b) Switching field versus the particle's diameter N 
for ip = 0, j = k s = 10 2 . 
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Fig. 10. Critical field versus the surface anisotropy constant for tp = 0, and different 
values of surface-to-core ratio of exchange couplings; N = 10. 

surface spins are aligned along their radial easy axes, and because of strong 
exchange coupling they also drive core spins in their switching process, which 
requires a very strong field to be completed. The value of k s where h c jumps 
up (e.g., = 1 for J s /J = 1) marks the passage from a regime where scaling 
with the SW results is possible (either with a -^-dependent or independent co- 
efficient) to the second regime where this scaling is no longer possible because 
of completely different switching processes. 

To estimate K s and the critical field, consider a 4 nm cobalt particle of fee 
crystal structure, for which the lattice spacing is a = 3.554/4, and there are 
4 cobalt atoms per unit cell. The (bulk) magneto-crystalline anisotropy is 
K c ~ 3 x 10~ 17 erg/spin or 2.7 x 10 6 erg/cm 3 , and the saturation magnetization 
is M s ~ 1422 emu/cm 3 . The switching field is given by H c = (2K C /M s )h c . For 
■0 = 0, kg — 1 and h c = 15, so H c ~ 6T. On the other hand, k s = 1 means 
that the effective exchange field experienced by a spin on the surface is of the 
order of the anisotropy field, i.e. zSJ/2 ~ 2K S . Then using J ~ 8mev we get 
K s ~ 5.22 x 10~ 14 erg/spin, or using the area per surface spin (approximately 
a 2 /8), K s ~ 5erg/cm 2 . For the case of ip — 7r/4, k c s ~ 0.2 and h c ~ 0.3, which 
leads to H c ~ 0.1 T and K s ~ 1.2 x 10~ 14 erg/spin or 1.2erg/cm 2 . 

4-2 Surface contribution to the energy of magnetic nanoparticles: NSA 

We calculate the contribution of the NSA [25] to the effective anisotropy of 
magnetic nanoparticles of spherical shape cut out of a simple cubic lattice. 
The effective anisotropy arises because of deviations of atomic magnetic mo- 
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Fig. 11. Magnetic structure in the plane z = of a spherical nanoparticle of linear 
size iV = 15 with L/J = 2 for the global magnetization directed along [110]. 

ments from collinearity and dependence of the energy on the orientation of the 
global magnetization with respect to crystallographic directions. We show [40] 
that the result is second order in the NSA constant, scales with the particle's 
volume, and has cubic symmetry with preferred directions [±1, ±1, ±1]. 

As shown many times before, as the size of the magnetic particle decreases, 
surface effects become more and more pronounced. In many cases surface 
atoms yield a contribution to the anisotropy energy that scales with particle's 
surface, i.e., to the effective volume anisotropy decreasing with the particle's 
linear size R as K v , e s = Kv + K s /R, as was observed in a number of ex- 
periments (see, e.g., Refs. [16], [41]). The 1/R surface contribution to K v>eS 
is in accord with the picture of all magnetic atoms tightly bound by the ex- 
change interaction while only surface atoms feel the surface anisotropy. This 
is definitely true for magnetic films where a huge surface contribution to the 
effective anisotropy has been observed. The same holds for cobalt nanoclusters 
of the form of truncated octahedrons [42] where contributions from different 
faces, edges, and apexes compete, resulting in a nonzero, although significantly 
reduced, surface contribution to Kv, e s- 

However, for symmetric particle's shapes such as cubes or spheres, the sym- 
metry leads to vanishing of this first-order contribution. In this case one has to 
take into account deviations from the collinearity of atomic spins that result 
from the competition of the surface anisotropy and exchange interaction J. 
The resulting magnetic structures (for the simplified radial SA model) can be 
found in [39], [43], [11] (see also Fig. 11 for the NSA). In the case L > J devi- 
ations from collinearity are very strong, and it is difficult, if not impossible, to 
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Fig. 12. Left: Effective anisotropy energy of the particle for different orientations 
of its global magnetization showing cubic symmetry of Eq. (41). Right: Differences 
of particle's energies between main orientations of the global magnetization vs the 
particle size in the scaled form for L/J = 0.1 and 0.01. The scaling is valid for 
N < J/L, and its violation for L/J = 0.1 is seen in the right part of the figure. 

characterize the particle by a global magnetization suitable for the definition of 
the effective anisotropy. For L <C J the magnetic structure is nearly collinear 
with small deviations that can be computed perturbatively in L/J^ 1. The 
global magnetization vector m [this is the same as the vector denned in (6)] 
can be used to define the anisotropic energy of the whole particle. The key 
point is that deviations from collinearity, and thereby the particle's energy, 
depend on the orientation of m , even for a particle of a spherical shape, due 
to the crystal lattice. To illustrate this idea, we neglect the bulk anisotropy 
and the DDI in the Hamiltonian (1). For a sc lattice Eq. (3) reduces to 

-iy(NSA) _ V- -iy(NSA) WNSA) _ L ^ 2 

' L an ~ ' L a,n,i > ' L an,i ~ r> A i<x°ia) \ OZj ) 

i Z a=x,y,z 



where Zi a = 0, 1, 2 are the numbers of available nearest neighbors of the atom 
i along the axis a. One can see that the NSA is in general biaxial. For L > 
and z ia = < z%p = 1 < z.- n = 2 the a-axis is the easy axis and the 7-axis is 
the hard axis. If the local magnetic moments Sj are all directed along one of 
the crystallographic axes a, then the anisotropy fields = —dH,Ai/dSi are 
also directed along a and are thus collinear with Sj. Hence, at least for L ^ J, 
there are no deviations from collinearity if the global magnetization mo is 
directed along one of the crystallographic axes. For other orientations of mo, 
the vectors and H^j are not collinear, since then at least two components 
of Sj are non zero, and the transverse component of H^j with respect to Sj 
causes a slight canting of Sj and thereby a deviation from the collinearity of 
magnetizations on different sites. This adjustment of the magnetization to the 
surface anisotropy lowers the energy. As we shall see, this effect is strongest 
for the [±1,±1,±1] orientations of mo. For both signs of L these are easy 
orientations, whereas [±1,0,0], [0, ±1,0], and [0,0, ±1] are hard orientations. 
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To solve the problem numerically, we have to fix the global magnetization of 
the particle in a desired direction z/ (\vo\ = 1) by using the energy function 
with a vector Lagrange multiplier A: 

F = H-N\-(y-v Q ), ^tP^T (33) 

1 &i I 



To minimize T we numerically integrate the evolution equations 
Sj = - [Sj x [si x Fi}\ , Fi = -dF/dsi 

A = dF/d\= -N [v - u Q ) . (34) 

starting from Sj = v = m and A= 0, until the stationary state is attained and 
an energy minimum is found. Our numerical results for the magnetic energy of 
spherical particles as a function of the orientation of the global magnetization 
are shown in Fig. 12 (left). Fig. 12 (right) shows differences between the basic 
directions [001], [Oil], and [111]. It is seen that AE/N ^const for N -> oo 
limit, i.e., AE scale with the particle's volume V ocAf. 

To analytically solve the problem in the continuous limit, we replace in Eq. (32) 
the number of nearest neighbors of a surface atom by its average value 

Zia ^z ia = 2- \n a \/ max{|n x |, \n y \, \n z \} . (35) 

Here n a is the a-component of the normal to the surface n. The surface-energy 
density can then be obtained by dropping the constant term and multiplying 
Eq. (32) by the surface atomic density /(n) = max{|n x |, \n y \, \n z \}: 

E s (m,n) = -— \n x \m 2 x + \n y \m 2 y + \n z \m 2 z . (36) 

At equilibrium the Landau-Lifshitz equation reads 

m x H eff = 0, H cff = H A + J Am. (37) 

For small deviations from collinearity one can seek for its solution in the form 
m(r) = m + ^(r,m ), ^ = |^| < 1 (38) 

where ip is the solution of the internal Neumann boundary problem 



A^ = 0, 

or 



= f(m,n) 

r=R 
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f = — 



J 



dE s (m, n) / dE s (m, n) 



dm 



dm 



m m 



(39) 



This equation can be solved with the help of the Green function G(r, r') (see 
Ref. [40]), and the final result for the second-order energy contribution is 



£2 = ^jj / d 2 rdVG(r,rO£ 5 (m,n)£ 5 (m,n'). 



(40) 



Taking into account the cubic symmetry and computing numerically a double 
surface integral one can write the result of Eq. (40) as 

S 2 k — - — (ml + m 4 + mf), k = 0.53465, (41) 
Jo 

where Jq = zJ = 6J. This defines the large-iV asymptotes in Fig. 12 (right) 
shown by the horizontal lines. 

The analytical results above is valid for particle sizes N in the range 

l«iV« J/L. (42) 



The lower boundary is the applicability condition of the continuous approxi- 
mation. Since the surface of a nanoparticle is made of atomic terraces sepa- 
rated by atomic steps, each terrace and each step with its own form of NSA 
[see Eq. (32)], the variation of the local NSA along the surface is very strong. 
Approximating this variation by a continuous function according to Eq. (35) 
requires pretty large particle sizes N. This is manifested by a slow convergence 
to the large- N results in Fig. 12 (right). The upper boundary in Eq. (42) is 
the applicability condition of the linear approximation in ip. For N > J/L 
deviations from the collinear state are strong, and the effective anisotropy of 
a magnetic nanoparticle cannot be introduced. 

As we have seen in Eq. (41), the contribution of the surface anisotropy to the 
overall anisotropy of a magnetic particle scales with its volume V oc N 3 ~ M . 
This surprising result is due to the penetration of perturbations from the 
surface deeply into the bulk. If a uniaxial bulk anisotropy K c was present 
in the system, perturbations from the surface would be screened at the bulk 
correlation length (or the domain-wall width) 5 ~ yj ' J/K c . Then for N > 5 the 
contribution of the surface anisotropy to the overall anisotropy would scale as 
the surface: £2 ~ (L 2 / J) N 2 5. For not too large particles, iV < S, contributions 
of both anisotropies to the anisotropic energy are additive and scale as the 
volume. If the bulk anisotropy is cubic, both contributions have the same 
cubic symmetry [c.f. Eqs. (4) and (41)], and the experiment should provide a 
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value of the effective cubic anisotropy different form the bulk value [42]. For 
the uniaxial bulk anisotropy, the two contributions have different functional 
forms. Even if the bulk anisotropy is dominant so that the energy minima are 
realized for m||e 2 , the surface anisotropy makes the energy dependent on the 
azimuthal angle ip. This modifies particle's energy barrier by creating saddle 
points and strongly influences the process of thermal activation [44]. 

For small deviations from the cubic or spherical shape, i.e., for weakly elliptic 
or weakly rectangular particles, there should emerge a corresponding weak 
first-order contribution E\ that would add up with our second-order contri- 
bution. For an ellipsoid with axes a and b = a(l + e),e 1, the anisotropy 
energy scales with the surface, £ 1 ~ LN" 2 ^em 2 z [cf. Eq. (41)], so that 

£ 2 LN 

i ~ 77 < 43 > 



can be large even for L/ J <C 1. 

The Neel constant L is in most cases poorly known. However, for metallic Co 
[45] quotes the value of surface anisotropy —1.5 x 10 8 erg/cm 3 , i.e., L ~ —10 
K. This is much smaller than J ~ 10 3 K, which makes our theory valid for 
particle sizes up to N ~ J/L ~ 100, according to Eq. (42). For this limiting 
size one has £2/^1 ~ l/ e that is large for nearly spherical particles, e C 1. 



4-3 Surface effects on the magnetization of a nanoparticle atT>0: MC 



In this section, we consider more realistic model of round-shaped (spherical or 
ellipsoidal) nanoparticle of simple cubic or spinel crystalline structure, uniaxial 
or cubic anisotropy in the core, and transverse or Neel anisotropy on the sur- 
face, described by the Hamiltonian defined by Eqs. (1), (2), (4), and (3). Using 
various techniques explained above, we compute the magnetization, induced 
and intrinsic, as a function of temperature and applied magnetic field, for dif- 
ferent values of the surface anisotropy constant and exchange coupling [4] [46]. 
We shall mainly focus on novel features stemming from the combination of 
anisotropy, field and temperature effects on the magnetization. 



4-3.1 Ferromagnetic particles with Neel's surface anisotropy 

In our simulations we consider the core ferromagnetic coupling J c and uniaxial 
anisotropy K c = 0.01J C . On the surface we adopt J s = J c , while the anisotropy 
is given by the Neel expression (3) with constant K s = 0.1J C . We ignore 
the DDI for simplicity. In particular, we are interested in how anisotropy 
affects the superparamagnetic relation (8) that has been shown to hold at 
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Fig. 13. Left: Core, surface and net magnetizations as functions of temperature. 
Right: Circles: Scaled graph for the induced magnetization m as a function of 
HNM/T, where M is the intrinsic magnetization. Full line: Langevin function 
L(x) = cothx — 1/x. 

all temperatures below T c for isotropic systems. As no analytical calculations 
are possible here, we resort to the Monte Carlo technique with global spin 
rotations. In Fig. 13 (left) we plot the core, surface and net magnetizations of 
a nanoparticle of 257 spins, as functions of temperature in zero magnetic field. 
These results do confirm what was obtained from the spherical model [see 
Fig. 1 (left)] for isotropic box-shaped systems, namely that boundary effects 
suppress the magnetization, and here we see that this effect is enhanced by 
the SA.A drastic effect of the SA is clearly seen in the field dependence of the 
scaled induced magnetization m as shown in Fig. 13 (right) at T = T^ IFA /8. 
In the presence of a strong SA the superparamagnetic relation (8) is no longer 
valid, and m/M strongly deviates from the Langevin function. 




Fig. 14. Left: Temperature dependence of the surface and core magnetizations for 
N = 3766, magnetization of the bulk system, and that of the cube with the spinel 
structure and periodic boundary conditions (PBC) with M = 40 3 . Right: Spatial 
variation of the net magnetization of a spherical nanoparticle of 3140 spins, as a 
function of the normalized particle radius, for r core = T/T^ OTe <C 1, and r core = 0.5, 

T core \- 
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4-3.2 Maghemite (^-Fe^Oj,) nanoparticles 

In this section we deal with the ferrimagnetic maghemite nanoparticles (7- 
Fe 2 3 ) having spinel crystalline structure, summarizing the results of [4], 
[46]. This time we include the DDI in the Hamiltonian of Eq. (1). The bulk 
anisotropy in such materials is cubic, still we are using a uniaxial anisotropy 
to simplify the study of effects that are of more interest to us here. We con- 
sider maghemite particles of various sizes (Af ~ 10 3 — 10 5 that correspond to a 
radius of 2-3.5 nm) and with the physical properties in the core (spinel crystal 
structure with lacuna, exchange and dipolar interactions, anisotropy constant, 
etc.) as those of the bulk, except that the anisotropy is taken as uniaxial. We 
use the TSA model with K s = 0.06 erg/cm 2 . All spins in the core and on the 
surface are identical but interact via different couplings depending on their 
locus in the lattice. We assume that the exchange interactions between the 
core and surface spins are the same as those inside the core. Although we 
treat only the crystallographically "ideal" surface, we do allow for a scatter 
in the exchange constants on the surface. In contrast, in Refs. [12] it was as- 
sumed that all exchange interactions are the same but there was postulated 
the existence of a fraction of missing bonds on the surface. 

In Fig. 14 (left), we see that the surface magnetization decreases more rapidly 
than the core magnetization as the temperature increases. Moreover, it is seen 
that even the (normalized) core magnetization per site does not reach its sat- 
uration value of 1 at very low temperatures, which shows that the magnetic 
order in the core is disturbed by the surface. This may also be due to lacuna in 
the spinel structure. In Fig. 14 (right) we plot the spatial evolution of the local 
magnetization from the center to the border of the particle, at different tem- 
peratures. At all temperatures it decreases with distance from the center. At 
high temperatures, the local magnetization exhibits a temperature-dependent 
jump, and then continues to decrease. This indicates that there is a radius 
within which the magnetization assumes relatively large values. This result 
agrees with that of [3] where this radius was called the magnetic radius. 



5 Conclusion 

We have demonstrated by different analytical and numerical methods the im- 
portance of accounting for the magnetization inhomogeneities in magnetic 
nanoparticles, especially in the presence of SA. The latter makes the magne- 
tization inhomogeneous even at T = and in general modifies the relation 
between the intrinsic and induced magnetizations. It also changes the mag- 
netization switching mechanism, since for large SA the particle's spins switch 
cluster-wise. For small SA we were able to calculate the spin canting in the 
particle analytically and to obtain a novel second order contribution to the 
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particle's overall anisotropy. It remains to generalize this result for nonzero 
bulk anisotropy. Another important task is to study dynamical implications 
of the many-body effects in magnetic nanoparticles. 
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